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ABSTRACT 



Context. The variable X-ray spectra of AGN systematically show steep power-law high states and hard-spectrum low states. The hard low 
state has previously been found to be a component with only weak variability. The origin of this component and the relative importance of 
effects such as absorption and relativistic blurring are currently not clear 

Aims. In a follow-up of previous principal components analysis, we aim to determine the relative importance of scattering and absorption 
effects on the time-varying X-ray spectrum of the narrow-line Seyfert 1 galaxy Mrk 766. 

Methods. Time-resolved spectroscopy, slicing XMM and Suzaku data down to 25 ks elements, is used to investigate whether absorption or 
scattering components dominate the spectral variations in Mrk 766. 

Results. Time-resolved spectroscopy confirms that spectral variability in Mrk 766 can be explained by either of two interpretations of principal 
components analysis. Detailed investigation confirm rapid changes in the relative strengths of scattered and direct emission or rapid changes in 
absorber covering fraction provide good explanations of most of the spectral variability. However, a strong correlation between the 6.97 keV 
absorption line and the primary continuum together with rapid opacity changes show that variations in a complex and multi-layered absorber, 
most likely a disk wind, are the dominant source of spectral variability in Mrk 766. 



Key words, accretion; galaxies: active 



1. Introduction 

Determination of the origin of X-ray spectra of active galac- 
tic nuclei has been mired in ambiguity, with different classes 
of candidate models predicting similar time-averaged spectra. 
A continuum in the form of a power-law is at the core of all 
popular models; a 'reflection spectrum' produced by scatter- 
ing from the surface of Compton-thick mat erial is also likely 
to be present, the putative accretion disk ( Guilbert & Ree s. 



igBStlLightman & Whitelll988l;lGeorge & Fabianlll991l) must 
exist near the continuum region and such reflection would 
produce a hard spectral component with strong Fe K emis 



sion, observable below 10 keV ( George & Fabianl 
' I 1 1 1 1 ' I ^^^^ 



2002t 



Tanaka et al. 



1991 



1995 



Zdziarski et al., 1995; Perola etal 
Nandra et al.. ,1997). Layers of Compton-thin gas are also 
thought to shroud the nuclear system, with strong and (rela- 



* Based on observations obtained with XMM-Newton, an ESA sci- 
ence mission with instruments and contributions directly funded by 
ESA Member States and NASA, also on observations obtained with 
Suzaku, a collaboration between ISAS/JAXA and NASA/GSFC, MIT 



lively) unambiguous evidence for the presence of zones cov- 
ering a large range of ionization-stat es and column densities 
(ICrenshaw. Kraemer & Georgel 1200 3h . If 'competing' models 
have the same constituents then what makes them different? 
The principal distinction is the degree to which reflection pro- 
cesses dominate versus absorption. Also of general interest is 
whether we can detect the signature of relativistic blurring in 
the reprocessed components. In practice, models dominated 
by blurred reflection are indistinguishable from those com- 
posed of layers of complex absorption. It has proved diffi- 
cult to progress with current data, where one typically has a 
short baseline for the observation and is limited to model-fitting 
time -averaged X-ray spectra with only mod est energy resolu- 



tion (iReeves et all 120041; iTurner et al.ll2005h . 



A possible way forward is to use the observed spectral 
changes over time to break the ambiguity. Rather than simply 
fitting models to time-averaged data we can attempt to disen- 
tangle components that arise on different physical scales and in 
different zones by their differing patterns of spectral variability. 
Marked spectral variability is indeed observed in local AGN 
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down to timescales of thousands of seconds, and the long ob- 
servations of these sources that have been performed recently 
yield good sampling of spectral variability and the possibility 
of isolating the dominant effects. 

The longer observations of AGN, that tend to only be per- 
formed once a mission is well-established, have turned out to 
be invaluable because exposure times > a day have pushed us 
into a regime where new results have been observable. For ex- 
ample, a long XMM observation of Mrk 766 has shown that 
the primary variable component of the spectrum has the form 
of an absorbed power-law, with accompanying ionized reflec- 
tion that varies with the continuum. The latter result first be- 
came evident with the discovery of the correlation between 
ionised line emissio n and continuum variations down to 10 ks 
( Miller et al. , 20061) . The se results also conc ur with the anal- 
ysis of the 2001 data bv iTurner et alj ( 120061) . where evidence 
for a periodic Doppler shift of ionised emission from a radius 
r ~ lOOr^ was found. 

These long datasets also inspired the first applications of 
principal compone nts analysis (PCA) to dat a from a single 
X-ray observation dVaughan & Fabianl |2004|) . By correlating 
flux at different spectral energies one can decompose data into 
variable constituent components, this method works well to 
break data into mathematically 'additive' components (emis- 
sion components, direct or reprocessed) but has a limitation 
because it cannot separate out the mathematically 'multiplica- 
tive' changes that would be caused , for example, b y varying 
absorption. The conclusion from the lVaughan & Fab ian (2004) 
PCA study of MCG-6-30-15 was that the spectral variabihty in 
that case was driven primarily by a power-law of varying nor- 
malization combined with a slowly-varying or constant hard 

component of reflection. 

2007h we inves- 



In the first paper in this series (iMiller et al 



tigated spectral variability in XMM observations of Mrk 766, 
a narrow-line Seyfert 1 g alaxy at redshift z - 0.0129 
dOsterbrock & Pog gel,ll985^ We appHed a PCA method based 
on singular value decomposition that preserved the full instru- 
mental spectral resolution. It was found that the source can 
be modeled adequately using just two additive spectral com- 
ponents: a variable component comprising a power-law and 
ionised, modestly-broadened Fe emission line together with 
relatively non-variable component(s) that shows a hard spec- 
tral form. 

The possible origins for the hard component have been dis- 
cussed at length in Paper I. A model where there is some scat- 
tered and self-absorbed fraction of the primary continuum and 
possibly a contribution from a distant neutral reflector can ex- 
plain the offset component shape and the spectral variability 
(hereafter the 'scattering' model). Alternatively, the hard flux 
could be a heavily-absorbed fraction of the continuum, and we 
denote this the 'partial-covering absorption' model ('p-cov ab- 
sorption'). In the latter model the PCA breakdown into additive 
components indicated that covering fraction changes produce 
the observed spectral variability, as the additive components do 
not change much in shape but rather in relative dominance as 
the source flux varies. Importantly, both these interpretations of 
the PCA work with out recourse t o rela tivistic blurring of any 
spectral component. iMiUer et alj (l2007b discussed fits that in- 



clude blurred reflection and concluded that significant blurring 
was not required to explain Mrk 766, although the source is 
consistent with the presence of such effects. 

Here, we present here a detailed analysis of time-resolved 
spectra of Mrk 766, fitting spectral models to individual time 
intervals. XMM has accumulated ~ 700 ks of data from this 
highly-variable Seyfert galaxy from 2000-2005. The XMM data 
discussed here have already been used to detect variability of 
the peak energy and flux of a somewhat broad, ionised Fe emis- 
sion line in the 6 - 7keV range (Miller et al., 2006). Here 
we show, for the first time, the PCA model fit to the data 
themselves, comparing the two physical interpretations of the 
PCA results presented by Miller et al. (2007), i.e. the scattering 
and p-cov absorption models, to understand whether scatter- 
ing/reflection or absorption effects are dominant. Of particular 
interest here is a follow-up of the time-dependent behavior of 
two absorption features found in the Fe K-band using PCA, a 
feature at 6.97 keV likely originates in H-like Fe while a feature 
at 7.3 keV is of uncertain origin. 

In this paper, we concentrate on the origin of spectral vari- 
ability in Mrk 766. In this sense, consideration of what we dub 
the 'scattering model' may be taken as including consideration 
of a relativistically blurred reflection model, since the scattered 



component may also be fitted by such a model (iMiller et al 



2007), although we note in section|6]how absorption line vari- 
ations may be used to distinguish relativistically-blurred from 
non-blurred models. Whatever the model fit to the 'hard' com- 
ponent, its temporal behavior is a key piece of information that 
may help us disentangle the possible explanations for its origin. 
In addition to this new treatment of the XMM data we present 
new Suzaku data allowing us to extend the analysis to ~ 50 
keV, and increase the temporal baseline for the study. 

2. The Data 

2.1. XMM-Newton data 

In this paper we utilise all XMM-Newton (hereafter XMM) 
European Photon Imaging Camera (EPIC) pn CCD 



(IStriider et all 1200 Ih data available for Mrk 766, cover- 



ing 2000 May 20 (science observation ID 0096020101), 2001 
May 20-21 (science observation ID 0109141301) and 2005 
May 23 (science observation IDs in the range 0304030[1 - 
7]01). The data reduction is described bv IMiller et alj (120061) 
and as in that work, the combination of some non-optimal 
modes for the Metal Oxide Semi-Conductor CCDs, photon 
pileup and inferior signal-to-noise led us to use only the pn 
CCD data in this analysis. 

To perform the time-resolved spectroscopy we found that 
even time-selections of 25 ks duration gave adequate signal- 
to-noise in all spectra (this yielded slightly lower effective ex- 
posure times in some cases, due to data gaps). This sampling 
gave a good match to the variability timescales evident in the 
data. Where datasets did not divide exactly into 25 ks bins, 
the final bin of an 'obsid' (orbital segment of data) was al- 
lowed to be slightly shorter or longer than this. This yielded 28 
spectral slices across all available data; 2000 May yielded two 
spectral slices, 2001 yielded five slices and 2005 yielded 21 
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spectral slices. With time-sliced spectra, despite our attempts 
to optimize the selections, the signal-to-noise at high energies, 
E > 7.5 keV, was still found to be low. To optimize the spec- 
tral fitting we binned all spectra consistently with energy bins 
equal in width to the half width half maximum (HWHM) spec- 
tral resolution of the EPIC pn instrument. The HWHM varies 
with photon energy, so the bin widths chosen also vary contin- 
uously with energy, from ~ 35 eV at £ ~ 1 keV to ~ 80 eV at 
E ~ lOkeV. 

The energy range chosen for the analysis initially is 1 .0 - 
9.8 keV, this focuses the study on gas closest to the nucleus 
and avoids modeling in detail emission/absorption lines below 
1 keV, some of which likely arise from gas relatively distant 
from the central engine (e.g. lOgle et alll2003l) . The truncation 
at 9.8 keV is where signal-to-noise reaches a very low level. 

2.2. Suzaku data 

Mkn 766 was observed by Suzaku from 2006 November 16 
00:34:54 to 2006 November 18 08:44:01, in XIS nominal 
mode. We used events files from version 1.3.2.6 of the Suzaku 
pipeline. 

The XIS data was reduced using tasks from v.6.1.2 of 
the 'ftools' package. It was screened with XSELECT to ex- 
clude data within the South Atlantic Anomaly (SAA) and data 
within 436 seconds of the SAA (using the T_SAA parame- 
ter). Additionally we excluded data with an Earth elevation an- 
gle less than 5° and Earth day-time elevation angles less than 
20°. Finally we selected good events with grades 0,2,3,4, and 
6 and removed hot and flickering pixels using the SISCLEAN 
script. This screening left us with efifective exposures of ~ 87 
ks for each of the XIS 0, 1, and 3 chips. Note that the spaced- 
row charge injection (SCI) was not active during this observa- 
tion. The XIS products were extracted from circular regions of 
2.9'while background spectra were extracted from a region of 
the same size off'set from the source (and avoiding the chip cor- 
ners with the calibration sources). The response and ancillary 
response files were then created using the tasks XISRMFGEN 
and XISSIMARFGEN, respectively. The source contributed 
98.6%, 95.3%, and 98.4% of the total counts in from the XIS 
0, 1, and 3 products, respectively. The orbit of Suzaku leads to 
the exposure noted above being spread out over a long, 200 ks, 
baseline. To obtain time slices with similar exposure to the 
XMM spectra we split the data into four even time periods cov- 
ering a baseline of 50 ks each, for spectral analysis. 

The PIN background events file was provided by the HXD 
instrument team, used in conjunction with the source events 
file to create a good time interval applicable to both the source 
and background. The background events file was generated us- 
ing ten times the actual background count rate, so we increased 
the efifective exposure time of the background spectra by a fac- 
tor of 10. We found the deadtime correction factor using the 
HXDDTCOR task with the extracted source spectra and the 
unfiltered source events files. After deadtime correction, we 
were left with an effective exposure of ~ 87 ks. To take into 
account the cosmi c X-ray background ( Boldt & Leiterl 1987 : 
Gruberet allll999h . Xspec version 11.3.2 lArnaudl (Il996h was 



used to generate a spectrum from a CXB model, normalized 
to the Suzaku field of view, and combined with the PIN simu- 
lated background file using MATHPHA to create an total back- 
ground file. The source contributed 6.7% of the total counts. 



3. The Models 

In this paper, as we have many spectral components in the fits, 
we refer to the ionised reflection that is viewed directly (and 
appears from the variability of the Fe K emission line to cor- 
relate with the directly-viewed powerlaw) as the 'ionised re- 
flector'. We denote the hard/off'set component, as the 'hard' 
or 'scattered' component. It is not clear whether the scatter- 
ing gas is Compton-thick in that case and the use of the term 
'scattered' allows an easy distinction between this steady com- 
ponent and t he ionised reflector that is correlated with the con- 
tinuum flux. Miller et al] ('2007*) also considered whether mod- 
els where the spectral variability was dominated by absorption 
changes might explain the data, concluding that the hard com- 
ponent could be described as an absorbed fraction of the pri- 
mary continuum and that changes in covering fraction could 
then explain the source behavior. As described in the introduc- 
tion we term this the 'p-cov absorption' model. Both 'scatter- 
ing' and 'p-cov absorption' models are consistent with the PCA 
provided the spectral variability is dominated by covering frac- 
tion changes in the p-cov absorption model. 

The general approach taken was to utilize the fits to the 
direct and offs et components from the application of PCA to 
the XMM data dMiller et all l2007b . then extend the analysis to 
include the broader bandpass Suzaku data, as a test of the model 
on data not used to create it. 



iMiller et al.l(l2007b fitted the spectral elements derived from 
PCA, to investigate the source in terms of additive spectral 
components. Here we fit the data directly using time-resolved 
spectroscopy and this allows us to test how well the two PCA 
interpretations fit each time period for Mrk 766, and to see how 
the different model parameters vary with time. This analysis 
also aims to expand upon the PCA by finding an explanation 
for the residual spectral variability that was not well modeled 
there. As noted above then, we apply two interpretations of the 
PCA decomposition: first we fit the 'scattering model' and then 
the 'p-cov absorption model' that uses varying covering frac- 
tion to explain the source. The d erivation of these m odels from 
the PCA is discussed in detail by Miller et al. ( 20071) . and so the 
two models to be utilized are simply summarized below. Here 
we use the ^ form of ionization parameter where ^ - L/nr^, 
L is the 1-1000 Rydberg luminosity, n the gas density and r 
the absorber-source distance. We quote ^ in units of erg cms ' 
throughout. 

3.1. The Scattering Model 

The scattering model (Fig.[TJ is that found from fitting the PCA 
hard component using the scattered primary continuum (whose 
slope was determined from fitting the PCA variable compo- 
nent). The full spectral model then consists of: 
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1. A power-law of constant slope T - 2.38 and variable nor- 
malization. 

2. An ionised reflector that varies with normalization tied to 
that of the power-law. 

3. A ~constant scattered component. This is represented by 
summing components 1 and 2 and then allowing a scalar to 
represent the scattered fraction. 

4. A complex layer of ionised absorption in front of the scat- 
tered component, represented by the cmbination of two lay- 
ers: N„(l) = 4.74 X 1022cin-2jog^(i) ^ -0.17, Nh(2) = 
3.97 X lO^-'cm"^, log^(l) - 3.05. This physical might arise 
if the reflecting/scattering surface has its own atmosphere, 
for example. 

5. A neutral reflector with no blurring using the xspec pexrav 
model with a Gaussian emission line representing Fe Ka at 
flie strength expec t ed for a slab of material inclined at 30° 
dOeorge & FabianL[l99ll) . 

6. A low-column, low-ionization absorber fully covering all 
components: Nh(3) = 2.77 x lO^'cm"^ log^(3) = 1.80 
covered all components during 2000-2001 while NhO) = 
2.41 X 10^'cm"^, log^(3) = 1.26 covered all components 
during 2005 and during the 2006 Suzaku observation. 

7. Two Gaussian absorption lines, allowed to float in normal- 
ization, with energies fixed at 6.97 keV and 7.3 keV (in the 
rest-frame of the host galaxy). 

Construction of this model assumed the highest state ob- 
served to date (in 2001) to represent the sum of directly viewed 
and reflected components; data from the lowest flux state ob- 
served in XMM data dominated the determination of the form 
of the hard component, although a small contribution from the 
direct component was allowed in that determination: the initial 
parameters values were chosen from the model fits to the PCA 
components. 

3.2. The P-Cov Absorption Model 

The p-cov absorption model (Fig. [Til is determined using sev- 
eral layers of ionized gas to explain the overall shape of the 
hard component, this complex absorption profile for the ab- 
sorbed fraction is consistent with many physical models. The 
model is composed of: 

1. A power-law of constant slope F - 2.18 and variable nor- 
malization 

2. An ionized reflector that varies with normalization tied to 
that of the power-law 

3. A large and complex layer of absorption composed of 
Nh(1) = 4.03 X 1022cm-2,log^(l) = 1.19 overlaying a 
second layer A?h(2) = 4.57 x lO^^cm"^ log ^(1) = 2.3 com- 
prised the complex absorber whose variable covering frac- 
tion accounts for much of the observed spectral variability. 
The column densities and ionisation parameters of these 
two layers were initially assumed to be constant through- 
out the time-resolved fitting for both models but later we 
found it necessary to allow Nh{1) to vary (Table 1). 

4. A narrow, neutral component of Fe Ka 

5. A low-column, low-ionisation absorber covers all compo- 
nents: Nh(3) = 5.38 x 102'cm-2,log^(3) = 2.27 for 2000- 




scattering gas 

l) + ^y2) 



absorbing gas N^3) 



partially covered direct scattered 

Fig. 1. A schematic representation of the inner regions of the 
models described in Section 3. The distant reflector is omitted 
from the scattering model, for clarity 

2001 A^h(3) = 3.91 X 102icm-2,log^(3) = 1.93 for 2005- 
2006. 

6. Two Gaussian absorption lines, allowed to float in normal- 
ization, with energies fixed at 6.97 keV and 7.3 keV (in the 
rest-frame of the host galaxy). 

Construction of this p-cov absorption model assumed the 
highest state observed to date to be uncovered with respect to 
this complex dense absorbing layer, although the low column, 
moderately ionized gas layer of component 5 covers the direct 
component (and all other components) at all times. The power- 
law index is flatter in this model than in the scattering model 
because in the latter the hard scattered component exists in all 
flux states, but it is absent in the p-cov absorption model. The 
high-state data thus determined the form of the underlying con- 
tinuum before the variable covering by absorption is applied. 
The fits then allow any additive combinations of the uncov- 
ered and covered fractions to fit each spectral slice. The fact 
that almost all shape parameters are fixed means the model has 
very few degrees of freedom and the method of construction 
is designed to isolate the key variables most efifectively. The 
uncovered and covered fractions are naturally anti-correlated 
with this construction, although they have not been strictly re- 
quired to always sum to the same total intrinsic flux (i.e. some 
intrinsic variation in continuum normalization is allowed but 
not required). 

3.3. Fitting the models 

Thus, following the results from the PCA, the slope of the 
power-law in each model was frozen at the best fit value found 
in each scenario. The ratio of power-law to ionized reflector 
was also frozen for each model and the normalization of the 
summed pair 1 plus 2 allowed to vary; this reprodu ces the 



known line/continuum correlation jMiller et al.L 120061) . It was 



found that the low column, low-ionization layer of absorption 
needed to have a different column density for the 2000-2001 
and the 2005 epochs, and the column density was fixed at the 
best-fit value for each of those two time periods, in each model 
run. 

The 2005 value derived from XMM was later found appli- 
cable to the Suzaku data and utilized in those fits. 
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Fig. 2. A p-cov absorption fit to a slice of the low-state XMM 
data (slice 9) illustrating the presence of the two absorption 
features, at 6.97 keV and 7.3 keV in the rest-frame of the host 
galaxy 

In both cases, the ionized reflector w as modeled using the 
'reflion' model of lRoss & FabianI ( l2005h . assuming solar abun- 
dances and that ^ was fixed at 1600 erg cms"', as determined 
from fitting the mean spectru m. Some blurring was allowed 
in the reflection spectrum via a iLaor et al.l (1199 ih convolution 
model, provided in xspec by the 'kdblur' function, this was 
a very modest blurring constrained to arise within 100-300 
from a disk inclined at 30° to the line-of-sight, and necessary 
to model the cr ~ 0.27 keV width of the ionized diskline. 

Layers of absorption were modeled using tables generated 
from xsTAR version 211n (IKallman et alll2004l) with solar abun- 
dances, and turbulent velocity 200 km s ' . The Gaussian model 
absorption lines at 6.97 and 7.3 keV were included in addition 
to the complex absorption components as these were not mod- 
eled by those absorption layers. The 6.97 keV line most likely 
arises from highly ionized gas and thus would have little or no 
signature at lower energies (unless the column density is enor- 
mous). The 7.3 keV line is of uncertain identification and thus 
it is most safely modeled in an isolated manner. 

For the scattering model, the initially free components 
were: the normalizations of the directly-viewed and scattered 
fractions of the summed powerlaw plus ionized reflector (i.e. 
1 plus 2). Additional free components were: the normalization 
of an additional neutral component of reflection, modeled us- 
ing the xsPEC pexrav model with a linked Fe Ka line of strength 
expected from a neutral reflector subtendin g 27r steradians to 
the i lluminating source and viewed at 30° dGeorge & Fabian , 
Il991i) and the normalizations of narrow Gaussian components 
to model absorption features at 6.97 and 7.3 keV (in the rest- 
frame of the host galaxy). 

For the p-cov absorption model the initial free components 
were: the normalizations of the unabsorbed and absorbed frac- 
tions of the spectrum (by design, these are anti-correlated in 
this model); and the Gaussian line normalizations, as before. In 
this model the neutral Fe line at 6.4 keV is not linked to a com- 
ponent of neutral reflection, but allowed to vary in isolation; 
this model assumes that line arises from transmission through 
the absorber and not from a neutral reflector. 

The absorption lines, while highly significant i n the mean 





XMM spectrum and in the PCA offset component (.Miller et al 



20071) are highly sensitive to the form of the spectral model, as 



Fig. 3. Scattering model (top) and p-cov absorption model (bot- 
tom) fits shown as green solid lines representing fits to two of 
the individual time slices of XMM-Newton data (slice 6, blue; 
slice 9, red). The agreement between data and model illus- 
trates the applicability of either model over the broad range 
of flux covered by the source during the XMM observations. 
The model lines are shown extrapolated down to 0.5 keV illus- 
trating how exclusion of the softest data from the fits has not 
compromised our ability to distinguish models; the soft data 
indicate a preference for the partial-covering absorption model 

they fall in a region of some complexity. For this reason, the 
line fluxes derived from each model differ slightly, and in the 
paper we present both sets of values for comparison. A close- 
up of the Fe K region is shown in Fig.|2]for slice 9 of the XMM 
data. 

4. Results 

4.1. XMM results 

Fitting first all 28 time slices of the XMM data over l-9.8keV in 
the observed frame, the scattering model gave a total goodness- 
of-fit - 4796 for 3752 degrees of freedom (d.o.f.) while 
the p-cov absorption model gave 3785 for 3752 d.o.f.. The fits 
confirm that despite having relatively few varying components, 
both our two physical representations of the PCA components 
do indeed describe the spectral variability to first order, and the 
models describe many spectral slices completely satisfactorily 
(Fig. [3). Successful first-order application of the two physical 
models to the data leaves us in a position to now determine the 
origin of the additional complexity that could not be accounted 
for from the PCA. This is a crucial step because we need to 
achieve a statistically superior set of fits before a meaningful 
comparison of models can be made. To this end, we examined 
the residuals to the worst fits, finding that allowing a change 
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in column density of one of the zones of absorption can ex- 
plain the remaining residuals. Allowing this extra degree of 
freedom significantly improves fits to both models. This makes 
sense, since the PCA could not handle multiplicative changes to 
the source spectrum (note that covering fraction changes alone 
would show up as additive changes in PCA, since the profile 
doesn't change but rather the amount of covered and uncov- 
ered source spectrum) and one might thus expect that a mul- 
tiplicative variation would be the origin of spectral variability 
not handled by PCA. The column density changes are a natural 
modification as both models incorporate significant columns of 
gas. Allowing the ionisation parameter of that zone to be free 
in addition to the column did not further improve the fits, so we 
did not allow ^ as a free parameter. The degeneracy between 
Nh and ^ means that while we have modeled these residual 
variations as column density variations as tabulated, it is im- 
possible to determine how much of this change is actually due 
to column changes and how much is due to ionisation changes. 

An absorber in the line-of-sight would respond to changes 
in the continuum flux with recombination timescales down to 
seconds for gas in the density regime likely for layers close to 
the nucleus, and even faster ionisation timescales. In the scat- 
tering model our absorption happens after scattering, so there 
may be significant time lags and we would not necessarily ex- 
pect to see ^ changes consistent with the continuum variation. 
In the p-cov absorption model the gas might be expected to re- 
spond to the continuum variability but it is not possible to un- 
ambiguously determine how much of the continuum variability 
is due to the covering variations and how much to intrinsic con- 
tinuum flux changes. Given these ambiguities we do not seek 
consistency between observed flux changes and the ionisation 
state of the absorbers. 

After freeing the column density for Nh(X) the fits were 
re-run with the column density of Nh{\) allowed to float in 
addition to the free parameters already noted. Freeing this col- 
umn improved both models: the scattering model was improved 
to an overall goodness-of-fit - 4274 for 3724 d.o.f. while 
the p-cov absorption model yielded 3488 for 3724 d.o.f. for 
the XMM data. The improvement to the fits in both cases in- 
dicates that column density changes are a significant compo- 
nent of spectral variability in Mrk 766 in either the scatter- 
ing or absorption scenarios. While the absorption model pro- 
vides a fit that appears statistically superior in the current model 
constructions, we note that the scattering model has been con- 
structed using basic model components to attempt to describe 
an interesting physical scenario. The recourse to this approxi- 
mation for the combined scattering and self-absorption effects 
that may occur in (for example) a disk wind, is due to the lack 
of a full self-consistent wind model at this time. For this rea- 
son, the simple fit statistic should not be a basis for ruling out 
this crudely-parameterized but interesting physical possibility. 

4.2. Suzaku results 

Extending our consideration to the Suzaku data we first cre- 
ated a light curve from the PIN data, and compared it to light 
curves in the 0.5-2 and 2-10 keV bands of XISO. Fig. |4] shows 
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Fig. 4. Top: XISO time series in 500 s bins in the 0.5-2.0 (red) 
and 2-10 keV (green) bands and from the PIN data covering 
~ 12-40 keV (black) Bottom: PIN light curve versus pre- 
dicted flux in the 12-40 keV band from the scattering (red) and 
absorption (green) models. 

those time series and it appears that while the XIS data vary 
by a factor ~ 6 from trough to peak, the PIN data are incon- 
sistent with a zero-lag coiTelation in variability. However, there 
are some uncertainties in the production of the PIN light curve 
that relate to the difficulty in producing an accurate subtraction 
of the instrumental background level: for this reason the PIN 
curve is not fitted to the model predicted rates, although these 
are illustrated in Fig.|4] 

When the PIN data are binned into the four time intervals 
chosen for spectral analysis the light curve is consistent with a 
constant flux and so we used the mean PIN spectrum to obtain 
a fit with each XIS time slice (Tables Al,2; Figs.|5]|6]l. 

Fitting the spectral data we found that three of four time 
slices were well explained in the same way as the XMM data 
(Tables A1,A2). However, the third time slice was a very poor 
fit to both models (reduced - 1.88 for the scattering model, 
X^- - 2.99 for the p-cov absorption model) until the second 
large gas column, Nh(2), was also allowed to float. Allowing a 
second layer of absorption to vary yielded significantly better 
fits to both models (Tables A1,A2) although yet more complex- 
ity is apparently indicated (Fig.|5]|6]l. 

Both physical interpretations of the PCA results provide 
reasonable fits to all XMM spectra. Tables Al and A2 show the 
fit parameters for each model, along with the goodness of fit for 
each slice. Fig.|7] shows the primary parameters for the fitted 
models, i.e. the flux of the two dominant spectral components. 
Comparison of Tables Al and A2 shows that the p-cov absorp- 
tion model provides the statistically superior fit to most time 
shces sampled from XMM while the scattering model is the 
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Fig. 5. Scattering model fits to the individual time slices of 
Suzaku data. Top: Red is slice 1, blue is slice 2; Bottom: Black 
is slice 3, Orange is slice 4. PIN data have been adjusted down- 
wards by the 13% flux correction determined from simultane- 
ous fitting of XIS and PIN spectra 



best fit to Suzaku data. However, the aim here is to determine 
whether scattering/reflection or absorption eff'ects dominate the 
spectral variability and, as noted above, if our template models 
are not perfect representations of the physical scenarios then 
we might be misled by a simple comparison of the statistic. 
To understand the nature of the source variations we need to 
consider which of the sets of results makes the most physical 
sense, and part of this involves consideration of whether any 
obvious parameter correlations emerge. 

4.3. Line Variability 

Fig.Eshows the absorption line fluxes plotted against time. The 
6.97 keV line appears to show significant variability while the 
7.3 keV line shows less indication of any variation; to quan- 
tify this we performed a test against a model representing 
a constant line flux. The fitted values of line flux were slightly 
diflerent for the two model applications because the line depths 
are critically dependent on the placement of the continuum. 

In the scattering model, the 6.97 keV line flux variations 
yielded - 66/3 1 d.o.f, variable at significance level ~ 2 x 
10""^, the 7.3 keV line flux variations gave - 28/31 d.o.f., 
consistent with a constant value. In the p-cov absorption model, 
the 6.97 keV line yielded - 71/31 d.o.f, also variable at 
significance level ~ 5 x 10"^, but the 7.3 keV line gave x^ - 
52/31 d.o.f., variable at significance level 0.01. In this case the 
evidence for variation in the 7.3 keV line turns out to be due 
to one deep line measurement during the low state in spectral 
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Fig. 6. P-cov absorption model fits to some of the individual 
time slices of Suzaku data, color coding and PIN data as for 
Fig.E] 

slice 9, exclusion of that single value shows the rest of the line 
fluxes to be consistent with a constant value. 

In spite of the dependence of the lines on the continuum 
model, the absorption line at 6.97 keV shows significant vari- 
ability in flux for both of the different parameterizations of the 
data. This strengthens our confidence in the robustness of the 
line variability result in this case. However the discrepant con- 
clusions on the 7.3 keV absorption line in the two scenarios 
leads us to conclude we cannot be sure of the variability behav- 
ior of this line component. However, the equivalent width of the 
7.3 keV line does vary systematically with the total source flux 
and is not consistent with being a line of constant equivalent 
width. 

The emission line at 6.4 keV also has different implied be- 
havior in the scattering and p-cov absorption models. In the 
former it is tied to the neutral reflector and measures of the 
line flux yield - 70/31 d.o.f (variable at significance level 
8 X 10"^) while in the p-cov absorption model it is not linked to 
any other continuum component and yields = 8.3/31 d.o.f, 
consistent with a constant value. The question of variability of 
this line component remains ambiguous, while the Fe emission 
varies in this regime it is impossible to disentangle the neutral 
line from the wings of the emission line from the ionized re- 
flector and the origin of the neutral fine cannot be conclusively 
determined. 

5. Parameter Correlations 

The presence or lack of compelling correlations between pa- 
rameters could help us distinguish between the scattering and 
partial-covering models. The achievement of reasonable fits 
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Fig. 7. Top: Direct (black, solid symbols) and scattered (red, 
open symbols) component amplitudes of the source flux for 
the scattering model. Bottom: Direct (black, solid symbols) 
and absorbed (red, open symbols) component amplitudes of the 
source flux from the partial covering model. Solid vertical lines 
show where long breaks in time have been omitted for clarity. 



(Tables Al, A2) means reliable parameters can now be ex- 
tracted and we can investigate the correlations between them 
for each of the models. 

First we plot the most robustly-measured parameters, the 
directly-viewed and scattered or absorbed fluxes. Fig. [T] shows 
the directly viewed flux (powerlaw plus correlated ionized re- 
flector) for both models (black line) compared to (red) the 
scattered fraction (top panel) and absorbed fraction (bottom 
panel). The flux history of the source reveals how the 2000- 
2001 epochs were dominated by a direct view of the power- 
law continuum and associated ionized reflector. The 2005 ob- 
servation found Mrk 766 at a very low flux rising later in the 
observation. In 2006 we again found the source in a fairly low 
state but with a rapid rise evident towards the end of the Suzaku 
observation. 




time /ks 

Fig. 8. Absorption line fluxes for lines at 7.3 keV (red, open 
symbols) and 6.97 keV (black, solid symbols) from the scatter- 
ing (top) and partial-covering (bottom) fits. 



Overall the steady level of the scattered flux in Fig.|6]is con- 
sistent with the PCA, that analysis indicated the scattered com- 
ponent to vary within bounds +38% of the mean, although the 
component is higher than expected at the start of the 2000 ob- 
servations. In the partial-covering model, some anti-correlation 
of directly-viewed and absorbed fractions is a consequence of 
the model construction. 

5.1. Absorption Line Correlations 

Next we investigate the relationship between the absorbed 
or scattered fraction, and the narrow absorption lines. The 
strength of the correlations between modeled quantities is as- 
sessed using Spearmann's rank correlation coefficient, assum- 
ing data points to be statistically independent, along with coiTe- 
lation parameters for all other combinations of variables in both 
models. The red-noise power-spectrum (see lMiller et al ] l2007h 
of the intrinsic source variations causes some degree of corre- 
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Table 1. One-tailed significance level of correlations between 
parameters in the scattering model (upper right quadrant) and 
p-cov absorption model (lower-left quadrant) fitted to XMM 
and Suzaku data and after allowing column variability for 
NhW- Anticorrelations are indicated by minus signs in front 
of each value of significance level. 
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Statistical uncertainties on the line flux values are sufficiently 
large that this red-noise effect has only a limited influence, al- 
though we should bear in mind that the effect may be present at 
some level and we should use the significance values as being a 
relative measure of which relationships are the most important 
rather than being an absolute measurement of significance. 

Considering first the scattering model results (Table[Tl up- 
per right quadrant), the flux of the 6.97 keV absorption line ap- 
pears to be anti- correlated with the directly-viewed component 
and not correlated with the quasi-constant scattered component 
(Figs.|9l[l0ll. 

In the p-cov absorption model parameterization the 6.97 
keV absorption line is anti-correlated with the direct flux or 
correlated with the absorbed flux, with a slight preference for 
the former scenario (Figs.l9l[T0]i. As the uncovered and covered 
fractions are anti-correlated in this model, it is hard to discern 
which has the primary correlation with the absorption line. 

The 7.3 keV absorption line flux does not show a corre- 
lation with anything in either model (Figs.l9l[T0li. However, 
the line equivalent width is not consistent with being constant 
when measured against the varying spectral component, but is 
consistent with a constant value against the scattered or ab- 
sorbed components in the two models. 



I typical error ^ 
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Fig. 9. Absorption line fluxes versus the hard component in the 
scattering model (top two panels) and p-cov model (bottom two 
panels) for XMM (dots) and Suzaku (triangles) data. 



lation between the modeled quantities, so the significance of 
correlations may be overestimated by this approach, but with- 
out more information on the power-spectrum of variations in 
each quantity it is hard to correct for this effect. In fact, the 



6. Discussion 

The ambiguity in interpretation of the hard X-ray spectra ob- 
served in the low flux state for many Seyfert galaxies has hin- 
dered progress understanding these sources. We have attempted 
to move forward by using a large dataset accumulated by XMM 
over 2000-2005, combined with data from the broad bandpass 
covered by Suzaku during a 2006 observation, for the highly 
X-ray variable narrow-line Seyfert 1, Mrk 766. 

Two models have been directly compared in application to 
time-selected spectra from Mrk 766. The aim of the compari- 
son has been to distinguish whether large-scale spectral vari- 
ability is best explained by changes in relative contribution 
from powerlaw plus ionized-reflector compared to the scatter- 
ing component, or by variations in the covering fraction of the 
absorber. It is very important to determine the dominant effects 
as several models that fit the mean spectra of Seyfert galaxies 
comparably well have very different physical interpretations. 

The most robust measurements are those of the directly 
viewed and scattered or absorbed fractions of flux. The obvi- 
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ous interpretations outlined by iMiller et alj (l2007h are briefly 
reviewed here to set the scene for further discussion. In the 
scattering model the hard emission component is steady, ex- 
plained by scattering from an extended region such as a disk 
wind that could also absorb scattered photons. PCA limits the 
flux variations in this component to +38% of the mean and 
this relatively low level of variability implies a size of at least 
several light days for this gas but with an indication of some re- 
sponse of the inner region of the wind to changes in continuum 
flux. 

The most likely candidate for the scattering gas is the 
highly-ionized inner zone that produces the 6.97 keV absorp- 
tion line. The scattered fraction of emission is a significant frac- 
tion of the total and so this zone likely has a large covering 
fraction, whereas the lower ionization layers further out prob- 
ably have too small a covering to produce a significant flux of 
scattered X-rays. Opacity changes observed in the data could 
be explained by variations in such a wind that aff'ect the degree 
of self-absorption of the scattered ph otons. 



Seyfert galaxie s generally (e.g. IVaughan & Fabian 2004 
Pounds et alll2004l) show the behavior exhibited by Mrk 766, 



i.e. a hard spectral form when the primary continuum flux is 
faint; this tells us something about spectral variability in the 
context of absorption models. The simple scenario of an intrin- 
sically variable continuum with additional variability caused 
by covering fraction changes from passage of clouds across the 
sight-line can be ruled out. The covering fraction would have 
to know about the continuum flux for the observed phenomena 
to be so widely seen, else we would see some highly absorbed 
high-state spectra both in general and also specifically within 
the long duration Mrk 766 observations: the observed variabil- 
ity would not be a natural consequence of such a 'mixed vari- 
ability' picture. Two variants of the p-cov absorption model 
seem plausible however, either the continuum drives apparent 
absorption changes by ionizing some of the gas (see later) or 
the continuum is intrinsically constant and the observed vari- 
ability simply a consequence of the changes in covering frac- 
tion. The p-cov absorption model fits the XMM data very well 
but is a less satisfactory fit to the Suzaku data; as the absorption 
model was built upon a PCA deconstruction of the XMM data, 
changes in the absorbing gas between those observations and 
the Suzaku epoch would naturally result in a poor fit to the lat- 
ter. The more we free the column and ionization-states of the 
layers, the better the Suzaku data are fitted, so in conclusion, 
the statistical quality of fits to Suzaku spectra could be inter- 
preted as evidence for absorption changes on year-timescales: 
these must be expected, when opacity changes on tens of ks are 
also evident in the data. 

One compelling alternative interpretation of the partial- 
covering behavior in the absorption model may be that the re- 
gion producing the hard X-ray continuum might vary in extent. 
In this alternative picture, when the putative coronal region pro- 
ducing X-rays is largest the continuum could be visible above 
gas structures that might normally hide the nucleus: this would 
give apparent changes in covering fraction as the comptonizing 
region expands and contracts. This could naturally yield a con- 
stant flux of absorbed continuum emission that is obscured and 
that dominates above 10 keV. In this case the correlated varia- 



tions in ionize d reflection an d power law continuum known in 
Mrk 766 (iMiUer et all l2007h could be attiibuted to the reflec- 
tor seeing more continuum photons when the coronal region is 
at its largest. However, perhaps a more standard interpretation 
of the 'p-cov absorption' picture would be to suppose that the 
absorbing material is an inhomogeneous disk wind, in which 
clumps of obscuring material pass between a central contin- 
uum source and the observer. A wind would naturally have the 
most ionized zones closest to the nucleus, consistent with these 
data. 

Considering the second-order spectral variabi lity (not ad 



dressed in the principal components analysis of iMiller et al 



2007): we find opacity variations in one of the absorption lay- 
ers can account for most of the observed second-order effects. 
Degeneracy between column density and ionisation parameter 
means we cannot isolate the origin of the opacity changes with 
current data but again, the opacity variations support the gen- 
eral picture of an important disk wind component, driving the 
observed spectral variability. 

It is interesting to note that both models presented here are 
broadly consistent with the soft-band data down to 0.5 keV. 
The partial-covering absorption model provides a particularly 
good fit and naturally explains the soft-band spectral curvature 
without recourse to any separate "soft excess" component. As 
the partial-covering model was not fit down to 0.5 keV, then 
the agreement of soft-band data with the extrapolation of the 
partial-covering model to soft energies, supports the applica- 
bility of that model. 

One important diagnostic result is that the abs orption line a t 



6.97 keV, discovered in the mean spectrum ( Miller et al.Ll2007h 



shows significant variability, correlated with other aspects of 
source behavior. Interestingly the absorption line itself is not 
modeled within any of the zones of gas used to describe the 
overall spectral curvature but, as evident in Tables Al, A2, the 
opacity changes themselves are anti-correlated with the flux 
of the direct component. All zones of gas may be respond- 
ing to continuum changes, or continuum changes themselves 
may be an artifact of variable absorption. This possibility has 
been considered previously: the tight correlation between the 
ionized reflector and continuum flux might be a consequence 
of variable occultation of both if Compton-thick 'bricks' of 
variable covering fraction occult continuum and reprocessed 
emissio n from within. Th e bricks would be implied to exist 
~ IQrg (IMiller et al.Ll2007 ) which is close to where the occulted 
line is estimated to arise (based on the emission line width). 
By similar arguments, opacity changes observed on compara- 
ble timescales driven by motion of clouds across the sight-line 
would indicate clouds at the same radial location. If the opacity 
changes originate as ionisation-state variations then the loca- 
tion is currently unconstrained. 

The existence of a strong isolated absorption line at 6.97 
keV yields the lower limit log^ ^ 4, ie. it comes from gas 
that is the most highly-ionized of all layers evident in the 2- 
10 keV band. It is clear that the other layers of lower ^ gas 
may exist along the same sight-line but at larger radial dis- 
tances from the central source: the highest-^ zone could not 
exist outside the other gas layers as shielding would prevent it 
achieving log^ ^ 4. Similar models have been discussed for 
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sourc es such as e.g. NGC 3516 (ITurner et al.l l2005b and NGC 
4151 dKraemer et all l2005b and thus it seems that spectral vari- 
ations driven by a time-variable disk wind may be applicable 
across the Seyfert population. 

The 6.97 keV absorption line is so strong that it is not 
possible to model the zone from which it originates with- 
out invoking a significant turbulent velocity for the absorb- 
ing gas, again, consistent with its origin very close to the 
active nucleus. Low turbulent velocities (few hundred km/s) 
predict much weaker lines even for columns densities of few 



xlO cm" as the line saturates quickly and the equivalent 
width predicted when line broadening occurs can be much 
higher. One must invoke a turbulent velocity ^ 2000 km s ' 
to fit the observed line at its deepest with a column that is 
Compton-thin; for a turbulent velocity of 3000 km s ' we find a 
fit for the line in the mean spectrum using a column of gas with 
Nh ~ 5 X lO^^cm"^. Variability in the absorption line at 6.97 
keV supports the general conclusion that rapid absorption vari- 
ation is very important : the presence of such a li n e strengthens 
previous suggestion s ( Nandra & Poundsl 1994; Turner et al. 



20051; iReeves et al.l |2004 that very large column densities of 



highly-ionized gas exist close to the nucleus, dangerously close 
to the conditions that cause most confusion with Fe emission 
features. 

The absorption-line correlations found here seem to present 
two possibilities: in the p-cov absorption model: the 6.97 keV 
line flux variations could be due either to being constant EW 
absorption in front of the fluctuating absorbed part, or to being 
over-ionized at times of exposure to high powerlaw flux. Fig.|9] 
shows that the 6.97 keV line flux does indeed correlate posi- 
tively with the flux in the absorbed component, lending weight 
to the first of these possibilities. There is however a stronger 
anti-correlation with the direct component (Fig.[TOl Table[T]i, 
which appears to show a sudden transition in line flux for di- 
rect component flux /2-iokeV ^ 10"" ergcm"^ s"'. We discuss 
below the difficulty of obtaining such a rapid transition purely 
from changing ionization parameter 

Compared with the 6.97 keV line variations, the 7.3 keV 
line behavior is less clear due to its model-dependency. One 
thing that is clear is that the feature is persistent, remaining 
evident in Suzaku data during 2006 and that it is not an arti- 
fact of any instrument. However, the identification of this line 
is unknown, there are no known transitions expected of this 
strength at such a rest- energy withou t a sub stantial complex of 



other lines appearing ( Kallman et al.L 12004 ). It is possible that 



this feature represents an ionized Fe edge whose shape is dis- 
to rted due to unresolv ed features close in energy, or as noted 
bv lMiller et al.l (l2007h the feature could arise in a high-velocity 
outflow such as from Fe XXVI Lya in gas with a bulk outflow 
velocity ~ 13, 000km"'. There is a suggestion that the 7.3 keV 
line variations are correlated with those of the 6.97 keV line, 
however, there is no evidence for a causal relation between the 
7.3 keV line and the direct component in the p-cov absorption 
model, and only weak evidence for correlation with the direct 
component in the scattering model (TablefTll. This lack of cor- 
relation is, in itself interesting. The equivalent width of the 7.3 
keV line is inconsistent with a constant value when measured 
against the variable spectral component in either model, but 



is consistent with a constant value when measured against the 
hard/steady component in either model: this may be associated 
with a quasi-constant continuum component, leaving the equiv- 
alent width variations attributed to changes in the flux of the 
direct continuum component. 

In the scattering model, the variation in 6.97 keV line flux 
would need to be explained as arising from a change in ion- 
ization state of the absorbing gas as the source brightens over- 
all. It appears that the decrease in line flux is associated with 
quite a sharp transition in total source brightness which would 
not be expected given the broad range of ionization param- 
eter over which any given ionization state generally exists 



(IKallman et al.L l2004l) . We can test this explanation by gen- 



erating model absorbers with xstar assuming that the line is 
Fe XXVI 6.97 keV and assuming that the ionization parame- 
ter varies with the total source brightness. As discussed above, 
the apparent lack of Fe Ka 6.7 keV in absorption indicates high 
ionization, log ^ > 4. As also discussed above, in order to ob- 
tain a high equivalent width without an extreme column den- 
sity some velocity broadening is probably required, and we test 
models with velocity broadening cr = 3000 km s"', a value that 
allows reconciliation of the measured line strength with physi- 
cal models, while staying within the velocity constraint yielded 
from the data (cr < 6000km s"'). The results are shown in 
FigEU where the equivalent width has been measured against 
the total continuum flux: i.e. it has been assumed that the ab- 
sorbing gas lies in front of all the source emission compo- 
nents. The radial distance and density of the absorbing mate- 
rial is unknown, so we have investigated the cases in the range 
n = 10*^-10'^ cm"^ and choose to adjust the ionization param- 
eter so that the maximum possible variation in line equivalent 
width is produced, corresponding to log^ ~ 5. 

It appears that indeed the rapid variation in equivalent width 
cannot be well reproduced simply by changes in source bright- 
ness affecting the ionization state of the absorber. It may be, 
however, that the particular absorbing zone lies only in front of 
the scattered component of emission, but that nonetheless its 
ionization state is still determined by the total source bright- 
ness. In this case the typical equivalent width measured against 
the scattered component is higher, but shows less variation. 
Although the ionization-variation hypothesis still does not pro- 
duce good agreement with the observed variation, such a pic- 
ture probably cannot be ruled out. This deduction does have an 
important implication for the nature of the scattering compo- 
nent. If the 6.97 keV absorbing gas has to be associated with 
this hard spectral component only, then that component cannot 
be significantly relativistically blurred, otherwise the absorp- 
tion line would also be blurred. But if the absorbing gas lies 
outside any such blurred region, the rapid variation in EW with 
source flux would not take place. This consideration seems to 
rule out relativistically blurred models for the scattered com- 
ponent, unless an alternative blurred model can be found that 
reduces the significance of the 6.97 keV absorption line. 

7. Conclusions 

Detailed time-resolved X-ray spectroscopy has confirmed that 
the spectral shape of Mrk 766 is directly coupled to the source 
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brightness. Two model interpretations of principal components 
analysis fit the time-sliced data well, with changes in spectral 
shape then attributed to the relative strength of the continuum 
and a self-absorbed scattered component, or to changes in cov- 
ering fraction of an absorber. In addition to this dominant mode 
of variation, significant temporal variations in opacity down 
to timescales of tens of ks are also required to fit the time- 
sliced data, as previously i ndicated frorn resid uals in the princi- 
pal components analysis (' Miller et al.L l2007l) . The previously- 
reported absorption line at 6.97 keV is likely Fe XXVI Ka from 
gas with log ^ ^ 4 and this varies correlated with the primary 
continuum flux, tracing the behavior of a highly-ionized gas 
zone and extending our understanding of the parameter-space 
spanned by the complex, multi-layered absorber in this Seyfert 
galaxy. 

In conclusion then, Mrk 766 shows strong evidence for 
stratified layers of absorption, possibly decreasing in ioniza- 
tion with increasing radial location. The gas layers have rapidly 
variable covering factor, opacity and ionization state and are re- 
quired to explain the spectral variability in this source, whether 
or not one allows the lowest flux state to be dominated by 
some residual scattered continuum component. Some scattered 
and reflected emission must be present to explain the obser- 
vations, especially the rapidly variable emission features cor- 
relation with continuum flux, that show up in the high state. 
These results, combined with evidence in the literature for sig- 
nificant bulk outflow veloci ties in X-ray absorbers in general 
(ICrenshaw. Kraemer & Geor ge. 2003) points to disk winds as 
dominant in shaping Seyfert spectral variations in the X-ray 
band. 
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Appendix A: Model fit parameters 

Tables Al and A2 give the model fit parameters for the scat- 
tering and absorption models applied to the XMM and Suzaku 
data. The model components are described in the main body 
of the paper. For each table, the values given in each time 
slice (column one) are: the 2 - lOkeV fluxes of the total, di- 
rect and scattered or absorbed components (columns 2-4); 
the absorption-Une flux in each of the 6.97 keV and 7.3 keV 
(rest-frame) absorption features (columns 5,6); the emission- 
Une flux in the 6.4 keV emission line (column 7); the Nh col- 
umn density (column 8); and the reduced noting that the 
XMM and Suzaku fits have difl'ering degrees of freedom. 




'PL 2-10 keV 




'PL 2-10 keV 



erg cm ^ s ' 




'PL 2-10 keV 



erg cm ^ s ^ 




'PL 2-10 keV 



/10- 



erg cm s 



Fig. 10. Absorption line fluxes versus the direct/unabsorbed 
component for the scattering model (top two panels) and p-cov 
model (bottom two panels) for XMM (dots) and Suzaku (trian- 
gles) data. 
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1 


' 1 ' 


1 ' 1 

- 







-11.2 -11 -10.8 -10.6 



l°9lo(fTOT2-10keV Al"?! =nn ^ s ^) 

Fig. 11. Variation of equivalent width of the 6.97 keV line, 
measured against the total source continuum, with the total 
source flux (points with error bars: for XMM (dots) and Suzaku 
(traingles) data.). Values are calculated within the scattering 
model. Also shown are model lines for the variation expected 
from XSTAR models of the Fe XXVI 6.97 keV line arising 
from ionisation variation alone, for column densities with solar 
abundance of A^h = lO^'-s, lO^^, lO^^-^, lO^^, lO^^-^ cm l 
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Scattering Model 



Time 


Total 


Direct 


Scattered 


6.97 keV 


7.30 keV 


6.4 keV 






Slice 


2-10 keV flux /lO-i' 


erg cm~^s'^ 


Line photon flux / lO * 


-2 -1 

cm 












2000 XMM 










1 


1.406 


0.782+0.010 


0.510+0.040 


1.1+2.0 


0.5+1.9 


9.4+1.0 




1.00 


2 


1.611 


0.955±0.024 


0.604+0.050 


1.7+2.9 


0.3+2.9 


4.5+2.0 


0.4!°:! 


1.03 








2001 XMM 








3 


2.724 


2.377+0.006 


0.296+0.036 


0.9+2.9 


2.2+2.8 


45+3.0 




1.56 


4 


2.560 


2.098+0.007 


0.344+0.030 


0.0+2.7 


0.0+2.7 


9.1+2.0 


o O+0.7 
-0.7 


1.29 


5 


2.715 


2.244+0.007 


0.401+0.030 


4.2+2.7 


1.4+2.7 


6.4+2.0 


3 6+° *' 


1.28 


6 


2.858 


2.340+0.004 


0.402+0.030 


1.2+2.9 


1.0+2.8 


9.7+2.0 


3 6+° *' 

■-'•"-0.7 


1.33 


7 


1.835 


1.329+0.007 


0.433+0.026 


3.4+2.5 


0.0+2.7 


6.2+2.0 


9 n+0.4 
^•-'-0.3 


1.01 








2005 XMM 








8 


0.548 


0.216+0.003 


0.318+0.017 


7.0+1.3 


4.6+1.2 


2.3+2.0 


7 Q+0.3 


1.20 


9 


0.920 


0.322+0.003 


0.555+0.017 


9.3+1.6 


5.4+1.5 


5.0+2.0 


3 5+"-^ 

-^■-^-0.2 


1.37 


10 


0.735 


0.260+0.003 


0.431+0.018 


8.3+1.4 


0.0+1.6 


4.3+3.0 


5.8!°| 


0.87 


11 


0.566 


0.187+0.004 


0.368+0.027 


3.9+2.4 


1.5+2.4 


1.5+3.0 


q o+o.i 

-1.0 


1.16 


12 


0.777 


0.377+0.004 


0.355+0.018 


6.5+1.5 


4.8+1.4 


48+3.0 


4 7+0.3 

^■'-0.3 


1.12 


13 


1.404 


0.959+0.005 


0.349+0.021 


5.3+1.9 


2.6+1.8 


8.6+3.0 


9 4+0.4 


1.26 


14 


1.217 


0.861+0.004 


0.297+0.020 


2.8+1.8 


3.3+1.8 


5.5+3.0 


T 7+0.5 
-'■'-0.1 


1.20 


15 


1.156 


0.815+0.004 


0.295+0.023 


4.1+1.9 


2.8+1.8 


4.5+2.0 


4 7+0.5 


1.10 


16 


1.301 


0.944+0.005 


0.307+0.021 


0.0+1.9 


0.0+1.9 


4.1+2.0 




1.38 


17 


1.554 


1.213+0.005 


0.299+0.025 


2.8+2.0 


4.4+1.9 


41+3.0 


4 Q+0.6 


1.31 


18 


1.444 


0.996+0.005 


0.418+0.022 


0.4+2.1 


0.0+2.1 


2.0+3.0 


Q Q+0.4 
^■^^-0.4 


1.10 


19 


1.264 


0.902+0.004 


0.337+0.025 


0.0+2.2 


1.6+2.1 


2.2+2.0 


4 Q+0.5 


1.11 


20 


1.482 


1.104+0.004 


0.345+0.023 


3.0+2.0 


0.9+1.9 


3.1+3.0 


4 5+"-5 

^■-'-0.5 


1.15 


21 


1.732 


1.391+0.004 


0.307+0.026 


2.1+2.1 


1.7+2.0 


3.2+3.0 


IT 1+0.6 
-^•■'-0.6 


1.04 


22 


1.953 


1.358+0.004 


0.566+0.024 


1.1+2.4 


2.0+2.3 


2.7+2.0 


9 Q+0.2 
^•^-0.3 


1.20 


23 


1.795 


1.188+0.004 


0.543+0.028 


0.0+2.7 


4.7+2.4 


5.7+3.0 


7 4+0.3 
^■^-0.3 


1.18 


24 


1.769 


1.349+0.004 


0.428+0.026 


5.6+2.2 


3.6+2.1 


0.3+3.0 




1.06 


25 


1.606 


1.229+0.004 


0.366+0.024 


3.4+2.0 


1.3+2.0 


1.4+3.0 


4 4+0.5 


1.20 


26 


1.247 


0.925+0.004 


0.284+0.023 


1.22+1.9 


4.3+1.8 


3.6+2.0 


4 4+0.6 


0.96 


27 


1.030 


0.532+0.003 


0.399+0.030 


9.2+2.6 


2.2+2.7 


0.9+2.0 


^■■'-0.5 


0.82 


28 


1.323 


0.921+0.002 


0.365+0.020 


1.9+1.9 


0.0+1.9 


3.2+2.0 


4 1+0.4 
^■-'-0.4 


0.87 








2006 Suzaku 








29 


1.325 


0.944+0.009 


0.297+0.020 


2.2+2.4 


4.7+2.4 


6.3+1.7 


7 Q+0.1 


1.17 


30 


0.768 


0.441+0.004 


0.215+0.016 


4.5+1.8 


5.0+1.8 


8.4+1.3 


7 Q+0.9 
^■'-0.6 


0.99 


31' 


1.238 


0.707+0.009 


0.407+0.020 


4.1 + 1.8 


2.8+1.8 


1.0+1.3 


^■■^-0.3 


1.16 


32 


2.344 


1.888+0.020 


0.443+0.050 


0.0+1.8 


0.0+8.0 


1.5+8.8 


^■'-0.8 


0.71 



Table A.l. Values of fitted parameters for the 'scattering' model, in each time slice, as described in the text. Notes: (1) XMM 
fits have 133 degrees of freedom, Suzaku fits have 145 degrees of freedom. (2) SUce 31 required A'h(2) free, with value 3.4;|;j ^ x 
lO^^cm-^ 
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Partial Covering Model 


Time 


Total 


Direct 


Absorbed 


6.97 keV 


7.30 keV 


6.4 keV 






Slice 


2-10 keV flux /10-" 


erg cm"^*"' 


Line photon flux / lO * 


-2 -1 

cm 


/lO^^cm-z 












2000 XMM 










1 


1.422 


1.079±0.010 


0.348+0.040 


2.2+2.0 


3.9+1.8 


0.0+5.0 


"•"-0.0 


1.04 


2 


1.617 


1.446±0.024 


0.173+0.050 


0.9+2.9 


1.9+2.7 


0.0+5.0 


"•^-0.4 


1.13 










2001 XMM 








3 


2.797 


2.793±0.006 


0.000+0.036 


0.0+2.8 


3.4+2.5 


4.5+5.2 


4 n+100.0 


0.87 


4 


2.556 


2.507±0.007 


0.035+0.030 


4.1+2.7 


2.1+2.5 


7.8+5.1 


69+113 


0.87 


5 


2.715 


2.682±0.007 


0.035+0.030 


0.0+2.7 


0.0+2.5 


4.6+5.0 


0.28+,^-^ 


0.79 


6 


2.859 


2.794±0.004 


0.061+0.030 


2.3+2.9 


3.8+2.6 


7.2+5.0 


0.3611} 


0.83 


7 


1.843 


1.610±0.007 


0.238+0.026 


0.0+2.5 


3.4+2.5 


0.0+5.0 


o.36:g:« 


0.83 










2005 XMM 








8 


0.550 


0.253±0.003 


0.307+0.017 


6.4+1.3 


4.4+1.1 


2.4+2.9 


9.34!l] 


1.01 


9 


0.930 


0.366±0.003 


0.578+0.017 


11.0+1.6 


9.4+1.4 


3.5+3.5 


1 47+"-'' 

^•^'-0.3 


1.29 


10 


0.738 


0.309+0.003 


0.435+0.018 


6.4+1.4 


0.0+1.5 


3.0+3.2 


S QQ+0.0 
• -0.0 


0.80 


11 


0.572 


0.217±0.004 


0.361+0.027 


4.0+2.4 


1.7+2.2 


0.9+4.4 


12.28-^°^; 


1.07 


12 


0.779 


0.447+0.004 


0.339+0.018 


5.8+1.5 


5.1 + 1.4 


4.8+3.3 


3.85+;;-^ 


1.11 


13 


1.405 


1.129+0.005 


0.280+0.021 


5.7+1.9 


4.8+1.7 


6.2+4.0 


o-o5:iii 


1.26 


14 


1.218 


1.006±0.004 


0.213+0.020 


2.5+1.8 


4.5+1.6 


4.8+3.7 


0.65!J:^ 


1.08 


15 


1.159 


0.955+0.004 


0.206+0.023 


3.7+1.9 


3.9+1.7 


5.2+4.0 




0.93 


16 


1.303 


1.097+0.005 


0.202+0.021 


0.0+1.9 


1.3+1.7 


3.0+3.7 


16*^-' 

"•^"-0.2 


1.01 


17 


1.552 


1.434+0.005 


0.122+0.025 


1.1+2.0 


5.3+1.8 


3.6+3.7 


3.50+1^ 


0.97 


18 


1.450 


1.157+0.005 


0.289+0.022 


0.3+2.1 


0.0+2.0 


0.0+4.0 


0.24il 


0.77 


19 


1.266 


1.063+0.004 


0.202+0.025 


0.0+2.2 


2.6+1.9 


1.2+5.0 


3.261\l 


0.93 


20 


1.477 


1.293+0.004 


0.187+0.023 


1.9+2.0 


2.4+1.8 


1.2+4.0 


1 36+^-'^ 


0.93 


21 


1.731 


1.641+0.004 


0.092+0.026 


0.1+2.1 


3.0+1.9 


1.4+3.8 


4 85+'-' 


0.79 


22 


1.960 


1.579+0.004 


0.386+0.024 


2.1+2.4 


5.1+2.1 


1.2+4.0 


"•"-^-0.0 


1.06 


23 


1.797 


1.391+0.004 


0.410+0.028 


0.5+2.7 


7.8+2.3 


2.9+4.0 


"•"-^-0.0 


1.12 


24 


1.774 


1.590+0.004 


0.193+0.026 


4.0+2.1 


4.5+2.0 


0.0+5.0 


3 16*^-^ 

J-i"_2.o 


0.78 


25 


1.609 


1.434+0.004 


0.180+0.024 


2.4+2.0 


2.9+1.9 


0.0+5.0 


"•"-^-0.0 


0.94 


26 


1.249 


1.089+0.004 


0.162+0.023 


0.5+1.9 


5.3+1.7 


3.1+5.0 




0.69 


27 


1.031 


0.626+0.003 


0.408+0.030 


8.9+2.6 


3.7+2.5 


8.3+5.0 


2.82+551 


0.82 


28 


1.325 


1.073+0.002 


0.252+0.020 


1.5+1.9 


1.3+1.8 


1.6+3.6 


0.83!^:° 


0.74 








2006 Suzaku 








29 


1.316 


1.160+0.009 


0.167+0.020 


4.5+2.5 


7.9+2.2 


2.5+5.2 


OKI 


1.38 


30 


0.764 


0.564+0.004 


0.215+0.016 


7.9+1.9 


7.9+1.7 


2.8+4.1 


0+"' 

"•"-0.0 


1.40 


31^ 


1.238 


0.801+0.009 


0.444+0.020 


5.8+2.3 


4.7+2.1 


0.0+4.9 


5 0+"* 

■^•"-0.7 


1.14 


32 


2.339 


2.264+0.020 


0.077+0.050 


0.4+8.1 


1.4+8.0 


0.0+14.0 


0+" '' 

"•"-0.0 


0.69 



Table A.2. Values of fitted parameters for the 'absorption' model, in each time slice, as described in the text. Notes: (1) XMM 
fits have 133 degrees of freedom, Suzaku fits have 145 degrees of freedom. (2) Slice 31 required A'h(2) free, with value 0.0 < 
0.5 X lO^^cm-l 



